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The influence of the tip-substrate bias induced electric field 
in a scanning tunneling spectroscopy experiment on colloidal 
InAs nanocrystals has been studied. Calculating the Stark 
induced splitting of the degenerate IPe state perturbatively 
within a particle-in-a-sphere model, revealed a possible ex- 
planation of recently publislapd experimental wave function 
mapping data by Millo et al.tJ. 



I. INTRODUCTION 

Mapping of the electronic wave functions in semicon- 
ductor quantum dots (QD) has recently become possi- 
ble using various experimental techniquescl u. Knowing 
the actual shape of the electronic densities contributes 
to a better understanding of the QD electronic structure. 
This knowledge is crucial with respect to the possible im- 
portance of semiconductor QDs as the ultimate building 
blocks of optoelectronic and nanoelectronic devices. 

Most recent, scanning tunneling microscopy (STM) 
measurementad also allow a wave function mapping of 
colloidal nanocrystals. The used InAs/ZnSe core/shell 
structures are spherical in shape resulting in atomic- 
like symmetries and degeneracies of the conduction band 
(CB) statetD'El. This leads, for example, to a six fold de- 
generate first excited state (here after referred to as the 
IPe state, where the subscription e denotes an electron 
rather than a hole state) such that a mapping of the IPe 
wave function should reveal the spherical superposition 
of all degenerate states. In ceptrast to this simple pre- 
diction the experimental datald show a torus like p,^2j^y2 
symmetry. 

We present a calculation of the CB states within a 
particle-in-a-sphere model taking into account the elec- 
tric field due to the applied STM voltage. The result- 
ing quantum-confined St§rk effect was studied earlier in 
three dimensional QDsat3. In this work we concentrate 
on the Stark effect induced degeneracy lifting of the first 
excited state which provides an explanation for the ob- 
served p,j.2j^y2 wave function symmetryEI at appropriate 
STM vohages. 



II. MODEL 

A sketch of the experimental setup in a scanning 
tunneling spectroscopy (STS) experiment on colloidal 



nanocrystal quantum dots is shown in Fig. ^ To obtain 
a tunnel spectrum the tip is positioned above a single 
nanocrystal attached to the substrate via hexane dithiol 
molecules. Keeping the tip-crystal distance constant the 
differential conductance as a function of applied voltage 
shows sharp peak^aO. For a detailed understanding of 
the experimental data it is thus crucial to know the dis- 
crete energetic spectrum of the QD since the obtained 



peak position 
dot structure 



Upeak are directly related to the electronic 



eUpeakiN, N+l)= j[E{N + 1, /i) - E{N, ly)] (1) 
E{N,fi = {n,}) = ^n,e, + V^°\{n,}) (2) 

i 

where the pre-factor 7 depends on the capacitive elec- 
trostatic geometry. For a very asymmetric tip-dot dot- 
substrate capacity distribution 7 is close to oneQ. The 
total energy E{N, {rii}) of a dot with N = J^i elec- 
trons, where denotes the occupation number of state i, 
is written as a sum of occupied single particle levels with 
energies and the total charging energy Vg"*. Hereby 
Vq°* includes both, the direct Coulomb interaction and 
all correlations. 




FIG. 1. Scanning tunneling spectroscopy of a single InAs 
nanocrystal. The InAs nanocrystal with a typical radius of 
a few nanometers are linked to a gold substrate by hexane 
dithiol molecules (DT). Trioctylphosphin (TOP) molecules 
form a ligand shell around the nanocrystal. At 4:.2K the tun- 
nel current is measured as a function of the applied voltage 
U between tip and substrate. 
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TS setup of Fig. |l| with voltages up to 
applied on a tip-substrate distance of a few 
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nanometers, the QD is exposed to a considerable elec- 
tric field. While a STM tip has to terminate in a single 
atom in order to achieve atomic resolution, the macro- 
scopic tip size is usually about one ordeun magnitude 
bigger than the here studied nanocrystalalJ. Other than 
a macroscopic metallic tip, a single terminating atom is 
not able to substantially focus the electric field. Over 
the nanocrystal size we can therefore assume the field 
between tip and substrate to be homogeneous in absence 
of the QD. Since the considered InAs nanocrystals are 
surrounded by ligands with a quite different relative di- 
electric constant compared to InAs, we modeled the QD 
as a jacketed dielectric sphere. Extending the text book 
calculation of the electric potential 4>inAs insidje-a dielec- 
tric sphere placed in a homogeneous field £ho-n^ to such 
a structure leads in spherical coordinates to: 

(j)InAs{r,9) = 

^e2£hoinrcos0 

2e2 + £162 -f 462 + 2ei + 2{^f[eie2 + £2 - ei - e^] 

with ei = 15.150 and £2 = 2. being the relative dielec- 
tric constants of InAs and the ligands respectively. i?i 
is the InAs core radius and R2 the total radius including 
the ligand shell (see Fift.. |l|) . As in the case of a dielectric 
sphere without a shelltll the core potential 4>inAs is still 
the potential of a homogeneous field. The field £hom oc- 
curring in (|^) is not directly accessible. It is obtained by 
equating the voltage U with the potential drop between 
tip and substrate of the inhomogeneous field outside the 
QD. 

The single particle spectrum of conduction band (CB) 
states was calculated usingj_a particle-in-a-sphere model 
with a finite potential wellEJ. Such a spherical confine- 
ment leads to size quantization comparable to the bulk 
energy gap of InAs. Hence non-parabolicity effects of 
the CB have to be taken into accouajt. We use an energy 
dependent effective mass approachEJ with 

m*{E) ^m*{Q)[l + E/Eg] (4) 

where m*(0) is the bottom CB effective mass and Eg 
the bulk energy gap. T|b£ validity of this approach has 
been checked by Bryantll3 for a dot radius of 3.2 nm by 
comparison to a 8 band calculation including CB-valence 
band coupling. 

Using the obtained electric potential inside the QD our 
model Hamiltonian is: 

with V{r) being the potential step with step height cor- 
responding to the work function of InAs towards air. The 
ligand shell is not included in the potential V{r). 

For calculating the QD states in the field we treated 
the electric potential perturbatively. The energies are 
obtained using second order perturbation theorytll. 



III. RESULTS 

Without an electric field the Hamiltonian (||) separates 
in an angular and a radial part. The angular Schrodinger 
equation is solved by spherical harmonics, and the radial 
Schrodinger equation by spherical Bessel functions ji ip^ 
side the well and spherical Hankel functions hi outsidetS 
the well. The continuity conditions at the potential step 
lead to a set of transcendental equation determining the 
energy levels: 

ahi{il3Ri) [lji-i{aRi) - {I + l)ji+i(ai?i)] = 
if3jiiaRi) [Ihi^iiiPRi) -{1 + l)hi+i{iPRi)] (6) 

with a = V2m*E/h and (3 = y^2m*(y - E)/h. The 
numerical solution of the first two energy levels for a dot 
radius of 3.2 nm is shown in Tab. 0. Other than in 
hydrogen the allowed orbital quantum numbers are not 
restricted by the principal quantum number. Hence the 
first exited state has the quantum numbers n = 1 and I = 
1. Owing to the spherical harmonics this state is 3-fold 
degenerate in the three magnetic quantum numbers m = 
— 1,0 and 1 (see Fig. ^ and These unperturbed wave 
functions of (||) are now used to calculate the effect of the 
electric potential, (jiinAs using second order perturbation 
theory. 
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FIG. 2. Stark splitting of the IPe level in a 3.2 nm radius 
InAs nanocrystal as a function of the applied tip-substrate 
voltage. For calculating the QD potential the ligand shell 
was taken to be R2 — Ri = 0.5 nm thick whereas the tip-shell 
distance was taken as another d = 0.5 nm (see Fig. ^. 3eV 
was-.used for the depth of the confining potential as obtained 
bycj using a pseudopotential approach. To get the energg 
dependent mass, the InAs bulk energy-igap Eg = 0A2eVt3 
and the effective mass m*(0) = 0.0239li3 in units of the free 
electron mass, were used. 

Although the ISe state also shows a Stark effect we 
concentrate here on the first exited IPe state, with re- 
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gacd to the experimental results obtained by Millo et. 
al.cJ. In contrast to hydrogen, the first excited state does 
not show a linear Stark effect, owing to the lack of s-p de- 
generacy in such a spherical well. The IPg degeneracy is 
lifted by the quadratic Stark effect, such that the energy 
of the lPe{m ± 1) wave functions oriented perpendicular 
to the field (bottom row in Fig. |^) are lowered compared 
to the lPe{m — 0) wave function oriented along the field 
(top row in Fig. a). This behavior can be qualitatively 
understood by looking at the electronic densities. In an 
electric field energy is gained by moving the electronic 
densities in field direction. The lPe{m ± 1) densities 
shown in the bottom row of Fig. ^ can move quite freely 
in that direction with only a small increase in confining 
energy. Hence a strong dipole moment is induced lead- 
ing to a pronounced Stark effect. On the other hand the 
cost in potential energy for the lPe(m = 0) density by 
moving into the same direction is much higher due to the 
close confinement potential step. This leads to a smaller 
induced dipole moment as can be seen in the top row of 
Fig. H, and therefore to a smaller quadratic Stark effect. 
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FIG. 3. Electronic density of the envelope wave func- 
tions calculated within a particle-in-a-sphere model. The 
densities illustrate qualitatively the different Stark effect on 
the IPe state, although they were calculated in an infinite 
spherical well neglecting the non-parabolicity of the conduc- 
tion band. Top and bottom left show the IPe(m ± 1) and 
XP^im — 0) densities without the presence of an electric field 
respectively. Whereas top right and bottom right show the 
corresponding densities in a homogeneous field of 0.18 V/nm 

The quantitative tip-substrate voltage dependence of 
the IPe energy level is shown in Fig. ^ for a QD ra- 
dius of Ri = 3.2 nmpi At an applied voltage of 1.4 V 
as used in experimentid the energetic difference between 
lPe{m = 0) and lPe(m± 1) is about 15 meV. This Stark 
induced energy splitting is experimentally resolvable and 



allows at an appropriate voltage tunneling into the ener- 
getically lower lPe{m± 1) states without tunneling into 
the lPe{m — 0) state. The qualitatively different density 
distributions of those split states shown in Fig. ^ are ob- 
served in a STM experiment by Millo et. al.EJ. Therefore 
the obtained Stark induced degeneracy lifting of the first 
excited IPg state serves as a possible explanation for ex- 
perimentally observed mapping of the lPe{m ± 1) wave 
functions only (see discussion). 
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FIG. 4. The plots above show the electronic density of Fig. 
^viewed along the applied electric field. The XP^ira = 0) den- 
sity is shown left whereas on the right side the superposition 
of IPfiim = 1) and \Pc{m = — 1) is plotted. 



IV. DISCUSSION 

As shown in this publication the first IPg electron will 
occupy one of the energetically lower IPe (to ± I) states. 
The second electron is then going to occupy the other 
lPe(rn=F 1) state. This leads to a torus like electron den- 
sity seen in a wave function mapping experiment. Due to 
the exchange interaction the spins of both electrons line 
up. 

If all IPe states had been degenerate the third elec- 
tron would have followed Hund's rule by occupying the 
IPe (to — 0) state with same spin orientation as the first 
two electrons. This configuration shown on the left of 
Fig. H would have led to a spherical electronic density. 
But due to the Stark effect the IPe degeneracy is now 
lifted such that configuration A competes with configu- 
ration B shown on the right of Fig. ^. This configuration 
would lead to a torus like electronic density since no elec- 
tron is occupying the lPe(TO = 0) state. 

In experimentu a torus like electron density was found 
with three electrons occupying the IPg state, suggesting 
that configuration B has a smaller total energy. 

Configuration A has to pay extra energy for the third 
electron due to the energetically higher lPe{m — 0) level. 
On the other hand all spins are lined up such that this 
configuration takes advantage of the full exchange en- 
ergy. In configuration B the situation is just opposite. 
While the third electron occupies the energetically lower 
IPe (to — 1) level this configuration lacks the exchange 
energy of the third electron due to its reversed spin. 
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m: 1 -1 



IPe 

m: 1 -1 

Configuration A Configuration B 

FIG. 5. Two possible configurations for tliree electrons in 
the IPe states. 

To find out which configuration has the lower energy, 
we calculate the contribution of the first order Coulomb 
interaction to the total energy. Using the wave function 
of a infinite spherical well we estimate that the exchange 
energy of configuration B is approximately 13 meV higher 
than in configuration A. This is 2 meV smaller than the 
Stark shift of A and therefore configuration B is favored, 
in agreement with experiment. 

The first order Coulomb contribution for two electrons 
in the ISe state yields a charging energy q£ 56 meV which 
is about half of the experimental valuecl. One reason 
for this small value is the neglected interaction with the 
polarization charges at the QD's surface. These charges 
induced by the electrons inside the QD are expected to be 
large due to the high dielectric mismatch at the surface. 
While this interaction enhances the charging |eiiergy, it 
has only a small effect on the exchange energyEl. 

Although our calculations explain the experimental 
data, the close competion between Stark effect and 
Coulomb interaction requires a more thorough investi- 
gation. On the one hand the actual electronic structure 
of the InAs/ZnSe nanocrystal, and on the other hand the 
influence of correlations on the Coulomb energy in these 
few electron systems has to be taken into account. 

V. CONCLUSION 

We calculated the Stark effect perturbatively on the 
first excited CB state in an InAs nanocrystal within a 
particle-in-sphere model. An energy dependent effective 
mass was used to account for the non-parabolicity in the 
CB. The obtained degeneracy lifting of the first excited 
IPe state provides a possible explanation for experimen- 
tally observed mapping of lPe{m ± I) wave functions 
without lPe{m — 0) mixing. 
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TABLE I. This table shows the first two energy levels in an 
Ri — 3.2 nm InAs QD without an electric field. The energy 
difference between the ground a^d first excited states is in 
good agreement with experimentu. 



State 


Energy in meV 


E{m) - E{lSe) 




542 




IPe 


861 


319 
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